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Abstract 

We study the isospin asymmetry in the isoscalar (IS) excitations in the mirror nuclei 14 and C 
by using the Hartree-Fock(HF)+random phase approximation (RPA) linear response function the- 
ory with a Skyrme interaction to take into account the continuum effect properly. The asymmetry 
in the IS monopole, dipole responses is pointed out in the continuum near the particle threshold 
with respect to the excitation energy and the sum rule strength. On the other hand, no clear sign of 
the asymmetry is found in the giant resonance (GR) region. In the quadrupole case, the calculated 
strengths of the mirror nuclei show almost the same energy dependence from the threshold to the 
GR region. It is found that the transition densities of the monopole response show an extended 
halo structure near the threshold, while those of GR region show a typical radial dependence of 
the compressional collective mode without any halo effect. Contrary to the transition densities of 
the monopole response, those of quadrupole response do not show any sign of the extended fea- 
ture of wave functions neither near the threshold nor the GR energy region. Calculated strength 
distributions of the IS multipole states are compared with recent experimental data obtained by 
the multipole decomposition analysis of a inelastic scattering on 1 O. 

PACS numbers: 21.10.Hw,21.10.Sf,21.10.Pc,21.60.-n, 21.60.Jz, 27.20. +n 
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I. INTRODUCTION 



A fundamental feature of nuclear structure is associated with a basic symmetry between 
neutrons and protons in the nuclear interactions. The existence of a general symmetry 
between np, nn and pp interactions was first observed in the low— energy nucleon— nucleon 
scattering data. The charge independence of the interaction implies the conservation of the 
isospin symmetry in nuclei [1, 2]. The charge symmetry assumption, based on the equality 
of the nn and pp interactions , might give almost the same excitation spectra in nuclei of the 
isospin partners. The isospin symmetry, however, is violated by the electromagnetic interac- 
tion. Additional symmetry— breaking effects arise from the neutron— proton mass difference, 
and small charge symmetry and independence breaking forces in the strong interactions. 

The isospin impurity was studied by using microscopic models in refs. [3, 4] in relation 
with Gamow— Teller and Fermi (5 decays. It was pointed out that the isospin mixing increases 
rapidly as a function of proton number and becomes few % in the ground state of 100 Sn. 
However it will be less than 0.1% in light nuclei with the mass less than A=16 and thus the 
isospin can be considered to be a good and useful quantum number in the light nuclei. For 
example, in the A=13 mirror nuclei 13 N and 13 C with the isospin T=l/2, all states in the 
spectra show one to one correspondence up to the excitation energy Ex=10MeV within few 
hundreds keV accuracy. 

In heavier nuclei, it was expected that the isospin symmetry might be of little significance 
because of very strong Coulomb field. However, the experimental observation of the isobaric 
analog states (IAS) in medium heavy nuclei shows that the isospin is well preserved even 
under the influence of the strong Coulomb field [5]. In light nuclei, the examples of the 
large isospin asymmetry were found in the first excited J 7r =l/2 + states of A=13 nuclei, 
13 C and 13 N, and also in the J w =l/2 + states of A=17 nuclei 17 and 17 F. This is called 
Thomas— Ehrman effect [6] . This large asymmetry is considered mainly due to the reduction 
of the Coulomb energy of loosely— bound low—/ state with small centrifugal barrier. It should 
be noticed that a large anomaly was also found in the isoscalar (IS) magnetic moments in 
A=9 mirror nuclei experimentally [7] and discussed in relation with the Thomas— Ehrman 
effect [8]. 

The Coulomb energy might be reduced not only for the loosely bound states but also for 
the unbound resonance and continuum states. In this paper, we will study the continuum 
states in the T=l mirror nuclei 14 and 14 C in order to extract the isospin asymmetry in the 
spectra. The A=14 isospin partners are particularly interesting because of a large difference 
in the threshold energies of the mirror nuclei. Namely, the neutron threshold energy of 14 C 
is S n =8.18MeV, while the proton threshold energy in 14 is S p =4.63MeV. This large energy 
difference between the two nuclei may give rise to a substantial effect on the low multipole 
strength distributions with J 7r =0 + and 1~, especially near the threshold. The lowest 1~ 
state in 14 has been obtained much attention due to its importance on the astrophysical 
CNO cycle in the stars [9]. The main aim of this paper is to study the isospin asymmetry 
of the continuum response functions in the A=14 mirror nuclei both near the threshold 
and giant resonance region. To this end, we use the Hartree— Fock (HF) + random phase 
approximation (RPA) continuum response function with a Skyrme interaction to calculate 
the IS strength distributions of low multipole states in 14 and 14 C. The two-body spin- 
orbit and Coulomb residual interactions are not included in the RPA response. This paper 
is organized as follows. A brief review of the RPA response function is given in Section II. 
Calculated results are shown in Section III and compared with available experimental data. 
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A summary is given in Section IV. 



II. CONTINUUM STRENGTH IN 14 AND 14 C 



The RPA linear response theory is based on the time— dependent HF (TDHF) theory 
[10]. The RPA Green's function G RPA is obtained as a small amplitude approximation of 
TDHF and expressed as 

bp 

= (i_G(o)^)-i G (0). (1) 
bp 

where is the residual particle— hole interaction. The unperturbed Green's function 
is defined as 



j 6 occupied 



1 1 



uj + irj — h + 6i uj — irj + ho — Cj 



nviin (2) 



where ipi and e« are the eigenfunction and the eigenenergy of the HF hamiltonian ho- The 
operator equation in the r.h.s. of (2) is nothing but the one-body Green's function in the 
coordinate space representation. We can use the standard technique to solve the Green's 
function taking into account the coupling to the continuum [11]. 

The transition strength S(u) for the states above the threshold can be obtained from the 
imaginary part of G RPA as 

s(u) = 1< n i f x,T ^ i >i 2 5 ^ - 

n 

= If f Im{f x ' T \7)G K?A {r^-uu)f x > T (r)}drd? 

(3) 

where f x ' T (f) is the transition operator with multipolarity A and isospin r. Below the 
threshold, the strength is calculated from the residue of the poles in the real part of the 
response function. The transition operators are given by 

A I 

yA-o.r-o = S ^ l r 1~/= for isoscalar monopole strength, (4) 



a 

f X=2 ' T=0 = Y j r 2 Y 2 ^(r i ) for isoscalar quadrupole strength. (5) 

i=i 

The IS dipole operator X)i r i^i^(^i) exc it es the spurious state corresponding to the center- 
of-mass motion. Thus, we consider the next order term in the expansion of the spherical 
Bessel function j\=\{qr) for qr<^ 1, 

A 

f x=1 ' T=0 = ffYi^fj) for isoscalar compression dipole strength. (6) 

i=i 
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The modified IS dipole operator is also defined as 

A 



/, A=1 ' T=0 = E^ 3 -^) y M(n) (7) 



where r\ = 3 < r 2 > /5 subtracts the spurious component from the operator (6) [12, 13]. 
The energy weighted sum rule (EWSR) is defined to be 

= E i < n \f> > i 2 = \ E < > (8) 

n,fi fi 

where H is the hamiltonian for HF and RPA calculations. For the IS excitations, the Skyrme 
two-body interactions do not give any contributions for the EWSR [10]. Thus the EWSR 
can be expressed for the operators (4), (5) and (7) to be 

S x=0 = ^-A<r"> (9) 

= Ll A < ^ > < 10 > 

The transition density for an excited state, | n> , 

5p(f) = <ri |E^-n)|0> (12) 

i 

can be obtained from the RPA response, since the imaginary part of G RPA (f, r' ; E n ) near 
the resonance is proportional to Sp(r) Sp(r')* . The reduced transition probability is calcu- 
lated using the transition density (12) as 

B (\,r:0^n) = E l< n \ | >| 2 = El / Sp(f)f^(f)df\ 2 (13) 

n ft 

where we use an identity of the one-body transition operator 

The radial transition density, 5p\(r) , is defined by 

Sp(r) = Sp x (r) YZJf) . (14) 

III. DISCUSSIONS 

The Skyrme interaction SkM* is used to calculate HF potentials and RPA response func- 
tions. The HF calculations are performed with the filling approximation putting all nucleons 
in the HF orbits from the bottom of the mean field potentials. In 14 0, the proton number 
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Z=8 corresponds to the closed shell while the neutron number N=8 in 14 C is the closed shell. 
Moreover, the energy difference between the 1^3/2 and lpi/2 states for neutrons in 14 and 
protons in 14 C are more than 6MeV in the HF calculations. Thus, the pairing correlations 
might not play any important role in the RPA response in these nuclei [14]. 

The Skyrme interaction SkM* gives the neutron separation energy S Tt =8.95MeV in 14 C 
and the proton separation energy S p =5.70MeV in 14 0. These values are somewhat larger 
than the empirical ones S n (exp)=8.18MeV in 14 C and S p (exp)=4.63MeV in 14 0. In order to 
perform quantitative study of the threshold effect in the response, we enlarge the spin-orbit 
strength of the SkM* interaction to be W o =150MeV-fm 5 in 14 C and W / " =155MeV-fm 5 in 14 
from the original value V4 / o =130MeV-fm 5 . These modifications give S n =8.11MeV in 14 C and 
S p =4.66MeV in 14 which are very close to the experimental ones. 

The existing Skyrme parameters do not reproduce the empirical separation energies of 
the two nuclei 14 C and 14 with the accuracy which is good enough to study the threshold 
behavior of excitations in loosely-bound nuclei. It has been known that the particle- vibration 
coupling has a substantial effect on the single-particle energy , especially around the Fermi 
surface [15]. For more sophisticated calculations, one should make an effort to readjust the 
Skyrme parameters including the particle-vibration coupling. Since it is not the aim of the 
present study to readjust all the Skyrme parameters, we take a conventional approach to 
obtain the empirical separation energies of A=14 nuclei. 

The unperturbed and RPA isoscalar monopole responses of 14 C and 14 are given in Fig. 
1. There are two bound particle— hole excitations z/, 7r(lsi/ 2 — >2si/ 2 ) {y and it stand for 
neutron and proton, respectively, hereafter) at Ex~31MeV in 14 C, while the proton particle 
state of the 7r(lsi/2 —^Si/2) configuration is in the continuum in 14 and gives a sharp peak at 
Ex~27.3MeV in the unperturbed response. The discrete states of the unperturbed response 
, which are not shown in Fig. 1, couple with the continuum in the RPA response and 
appears as a sharp peak at Ex~30MeV in 14 C and also a peak at Ex~28.5MeV in 14 0. The 
unperturbed peak at Ex~27.3MeV in 14 is integrated in the giant resonance peak in the 
RPA response. 

The large threshold strength can be seen just above the proton threshold in 14 and also 
above the neutron threshold in 14 C. The energy of the threshold peak in 14 is more than 
one MeV lower that that in 14 C due to the lower proton separation energy S p . One can 
see that the RPA correlations do not make any appreciable changes on the unperturbed 
response around the threshold peaks as was pointed out in ref. [16]. The monopole strength 
below 14MeV in 14 exhausts 10.1% of the energy weighted sum rule value (EWSR), while 
the monopole strength below 14MeV in 14 C has only 8.1% of the EWSR value. The giant 
resonance peaks appear at Ex~22MeV both in 14 and 14 C where RPA correlations enhance 
the monopole strength. The integrated monopole strength up to GR region (Ex<25MeV) 
is similar in the two nuclei 49.6% in 14 and 48.5% in 14 C, respectively. The transition 
densities of the threshold peaks are shown in Figs. 2 (a) for 14 C and (c) for 14 0. The 
dominant contribution comes from the neutron excitation to the continuum in 14 C, while 
the proton continuum excitation has the major contribution in 14 0. The neutron density 
shows a halo-type long tail extended more than lOfm in the continuum excitation of 14 C, 
while the proton density is rather compact and gives a minor contribution. The same halo 
effect appears in the proton transition density in 14 near the threshold, while the neutron 
one is very small and compact. Contrarily, the transition densities of GR region of 14 C and 
14 in Figs. 2 (b) and (d), respectively, have substantial contributions from both protons 
and neutrons. The radial dependence of the two contributions is very similar showing the 
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peak and the node at almost the same radii. The ratio between the proton and the neutron 
contributions around the peak at r~4fm is almost the same as that of the proton and the 
neutron numbers of each nucleus. This is one of the typical feature of the collective IS 
GR [17]. In general, one can see in Fig. 2 that the proton contribution to the transition 
densities in 14 is almost identical to the neutron one in 14 C both in the radial dependence 
and in the ratio to another contribution; the ratio of proton to neutron contributions in 
14 is very close to that of neutron to proton ones in 14 C. The IS monopole strength in 
14 was observed recently by the multipole decomposition analysis of a inelastic scattering 
experiments [18]. The IS peak is observed at Ex~7MeV and the integrated strength up to 
13.8MeV exhausts about 8% of the EWSR. These experimental data can be compared with 
the calculated peak at Ex~8MeV which exhausts 10% of the EWSR value. The observed 
strength in the GR region is about 38% of the EWSR up to Ex=25MeV in the a inelastic 
scattering. This value is close to the integrated strength of RPA response up to Ex=25MeV, 
i.e., 49.6% of the EWSR. 

The isoscalar compressional dipole responses of 14 C and 14 are given in Fig. 3. In 
order to subtract accurately the spurious component, I took a small mesh size 0.075fm and 
the reference radius 15fm. Then more than 95% of the spurious center of mass motion was 
eliminated from the physical states [19]. In Fig. 3, we adopted the dipole operator (7) where 
the spurious component is further subtracted from the operator. The dipole response of 14 
shows a sharp peak just above the threshold at Ex~6.1MeV and a broad bump peaked at 
Ex~9MeV. In 14 C, a discrete 1~ state is found at Ex=7.1MeV together with a broad peak 
at Ex~llMeV. In the two nuclei, substantial strengths are found also above Ex=15MeV. 

The observed lowest 1" state is at Ex=5.17MeV in 14 and at Ex=6.09MeV in 14 C, 
respectively [20]. It is remarkable that the calculated energy difference AE=lMeV be- 
tween the lowest 1~ RPA states in the two nuclei is as large as the experimental one 
AE(exp)=0.92MeV, although the observed 1~ states are found at somewhat lower exci- 
tation energies than the calculated ones. The strength of the threshold 1~ peak in 14 has 
2.5% of the EWSR value, while the discrete lowest 1~ state in 14 C has only 1.1% of the 
EWSR strength. The low energy peaks above the first excited states show also large differ- 
ences both in energy and in transition strength. Namely, the peak energy in 14 is 8.9MeV, 
while that is shifted to about 2MeV higher in 14 C with Ex(peak)~llMeV. The integrated 
strength up to Ex=12MeV is 17.3% of the EWSR in 14 which is much larger than that of 
the value in 14 C with 10.7% of the EWSR value. On the other hand, in the GR region, the 
two nuclei have almost the same transition strength, i.e., 23% of the EWSR between 12MeV 
< Ex<25MeV both in 14 and 14 C. 

The transition densities of several 1~ states are shown in Fig. 4; (a) and (d) for the 
lowest 1~ states, (b) and (e) for the threshold peaks, and (c) and (f) for the GR peaks. In 
order to subtract the spurious component from the transition density, we define a modified 
transition density 



where po is the total density of the ground state. The factor a is determined to satisfy the 
condition 



which is equivalent to evaluate the strength function with the operator (7) [13]. The tran- 
sition densities for the lowest V states show similar radial dependence in Figs. 4 (a) and 



op — a—— 

dr 
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(d) for 14 C and 14 Respectively, although the larger contribution comes from neutrons in 
14 C and from protons in 14 0. The transition density of the peak around Ex=11.0MeV in 
14 C (Fig. 4 (b)) is dominated by a large neutron contribution z^(lp 3 /2 — >2si/ 2 ) with a long 
tail extended until r=10fm, while that of the peak around Ex=9.0MeV in 14 (Fig. 4 (e)) 
with the dominant configuration 7r(lp 3 /2 — >Si/2) does not show any long tail. This difference 
might be caused by the Coulomb barrier for the dominant proton wave functions in 14 0. 
It is interesting to notice in Fig. 4 (f) that the proton and the neutron contributions have 
no node in the transition densities, but the IS one shows the node because of the different 
extensions of the radial wave functions. 

In ref. [18], the IS dipole strength in 14 is observed both in the low energy region and 
the high energy GR region. The experimental data show two peaks in the low energy region 
below Ex=10MeV. The integrated strength exhausts about 24% of the EWSR value below 
Ex=13.8MeV. The two peak structure and the similar sum rule value (18.5% of EWSR up 
to 13.8MeV) are found in the calculated results of Fig. 3. On the other hand, the observed 
dipole strength up to GR energy region (up to Ex=25MeV) amounts to be 70% of the EWSR 
value which is much larger than the calculated value of 40.5%. 

In the isoscalar quadrupole response of Figs. 5 (a) and (b), the discrete states are found in 
both nuclei at Ex=7.1MeV with B(IS)=94.6 fm 4 in 14 C and at Ex=7.2MeV with B(IS)=52.8 
fm 4 in 14 0, respectively. As far as the excitation energy is concerned, the isospin asymmetry 
is much smaller in the quadrupole state than the dipole states. This may be traced back to 
the difference of the particle configurations between the 1]" and 2+ states in the two nuclei. 
For the 2f states, the main configurations are z/(lp 3 / 2 — ^Pi/2) fc> r 14 and 7r(lp 3 / 2 -^lpi/2) 
for 14 C, respectively, and the lpi/2 states are bound in both nuclei. The bound particle 
states for for 2f excitations cause only a small energy difference between 14 and 14 C. On 
the other hand, the main configurations for the 1~[ states are 7r(lpi/ 2 — >Si/ 2 ) for 1 O and 
z/(lpi/2 — >2si/ 2 ) for 14 C, respectively. The particle state z/(2si/2) is bound in 14 C, while the 
state 7r(si/ 2 ) in 14 is the continuum state. This difference in the particle configurations is 
the origin of the large difference in the excitation energies of 1^ states between 14 and 14 C. 
The experimentally observed lowest 2 + state is at Ex=6.59MeV in 14 and at Ex=7.01MeV 
in 14 C. These data are consistent with the calculated smaller isospin asymmetry in the 
2 + states than that of the lowest l~ states. Namely, the observed energy differences are 
AE=0.92MeV for the 1" states, while that is AE=0.42MeV in the case of 2+ states. The 
calculated results show also the same trend. 

As is seen in Fig. 5, the threshold strength does not show any substantial enhancement 
in the quadrupole response, but the strong GR peaks are clearly seen around Ex=20MeV in 
the two nuclei. The EWSR of IS quadrupole response below Ex =25MeV exhausts 85.3% in 
14 C and 85.5% in 14 0. As far as the GR is concerned, the two nuclei show almost the same 
excitation spectra both in energy and in strength distribution. The transition densities of 
the lowest 2 + and GR states are shown in Fig. 6. It is interesting to notice that the proton 
amplitude is larger in the lowest 2 + state of 14 C, while the neutrons have the dominant 
contribution in 14 0. This is because the proton excitation (lp 3 /2 — ^Pi/2) is allowed in 14 C, 
while the neutron excitation (lp 3 /2 ^lpi/2) is allowed in 14 0. In Figs. 6 (b) and (d), 
the transition densities of quadrupole GR show a typical collective IS nature, i.e., the radial 
dependence is surface peaked, so called, Tassie-type and the neutron (proton) contribution is 
larger than the proton (neutron) one in the neutron-rich (proton-rich) nucleus 14 C ( 14 0). One 
can not see any substantial long-tail in the transition densities of the quadrupole response 
even in the low energy states. This is due to a larger centrifugal barrier for the quadrupole 
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response than the monopole and the dipole responses [21]. 
IV. SUMMARY 

We studied the isospin asymmetry in the IS multipole responses of A— 14 mirror nuclei 
14 C and 14 using the HF+RPA theory with Skyrme interaction to take into account the 
continuum coupling properly. We found that the IS strength distributions near the threshold 
show a large asymmetry between the two nuclei for the monopole and the dipole responses. 
On the other hand, the strength distributions in the GR region are very similar in both the 
peak energy and the sum rule value. It is seen a clear sign of the halo-type extension of 
the transition densities in the threshold monopole response. Contrarily, those of GR peaks 
do not show any sign of the halo effect, but give a typical IS collective features. For the 
quadrupole response, the two nuclei show very similar features not only in the GR region, 
but also in the lower energy region near the threshold. The calculated results of 14 O are 
compared with recent experimental data obtained by the multipole decomposition analysis 
of a inelastic scattering. The calculated monopole and the dipole strength distributions 
show good agreement with the empirical ones near the threshold. The empirical sum rule 
strength of the monopole states is also close to the calculated one in the GR region below 
Ex=25MeV, while the empirical dipole strength is reported to be larger than the calculated 
one in the high energy region between 12MeV< Ex <25MeV. In order to confirm further the 
isospin asymmetry of A— 14 nuclei, we need the experimental data of 14 C in the continuum 
region. 
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Figures 




FIG. 1: (Color online) Unperturbed and RPA IS monopole response S in 14 C and 14 calculated 
by the self-consistent response function theory with the Skyrme interaction SkM*. Unperturbed 
and RPA responses (3) are calculated by using the Green's functions (2) and (1), respectively, with 
the operator (4). See the text for details. 
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FIG. 2: (Color online) Transition densities (12) of the IS monopole states in 14 C and 14 0; (a) and 
(c) at the threshold peaks, and (b) and (d) at the GR peaks. The transition density is shown in 
an arbitrary unit. 
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FIG. 3: (Color online) IS compressional dipole response S in 14 C and 14 calculated by the self- 
consistent RPA response function theory with the Skyrme interaction SkM*. The dipole operator 
(7) is adopted in the calculations. The dashed— dotted line at Ex =7.75MeV shows a discrete 1~ 
state in 14 C. See the text for details. 
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FIG. 4: (Color online) Transition densities of the IS dipole states in 1 C and 14 0; (a) and (d) of 
the lowest 1~ states, (b) and (e) at the threshold peaks, and (c) and (f) at the GR peaks. The 
transition density is shown in an arbitrary unit. 



13 



140 
120 
100 



ffi 80 



■ (a) 


14 C IS quadrupole 




unpertubed A 




- 


— IS 






B(IS)=94.6fm 2 








1 '\ 
\ 1 ' 



10 15 20 

Ex(MeV) 



25 



30 



100 



80 



5 60 



40 



20 



(b) « 


IS quadrupole \ 




— unpertubed j 




— IS 




B(IS)=52.8fm 4 / 







10 15 20 25 30 
Ex (MeV) 



FIG. 5: (Color online) IS quadrupole response S in 14 C and O calculated by the self-consistent 
RPA response function theory with the Skyrme interaction SkM*. The operator (5) is used in the 
calculations. The solid lines show a discrete 2 + state in 14 C and a narrow resonance 2 + state in 
14 0. See the text for details. 
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FIG. 6: (Color online) Transition densities of the IS quadrupole states in 14 C and 14 0; (a) and (c) 
of the lowest 2 + states, and (b) and (d) at the GR peaks. The transition density is shown in an 
arbitrary unit. 
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